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ABSTRACT 

Momentum feedback from super-Eddington accretion offers a simple explanation for the observed M-a 
and M-M SD her relations bet ween supermassive black holes (SMBH) and the spheroids of their host galaxies. 
Recently Ferrarese et al. and Wehner & Harris observed analogous relations between the masses of central star 
clusters and their hosts. We show that stellar winds and supernovae from such nuclear clusters (NC) give 
similar feedback explanations for this case also, and we discuss the connection to the Faber-Jackson relation 
for the spheroids themselves. 

Subject headings: galaxies: star clusters, formation, nuclei 



1. INTRODUCTION 

It is well known that the masses of supermassive black holes 
(SMBHs) in the nuclei of early-type galaxies and late-type 
bulges correlate tightly with the velocity dis persions of the 
stellar spheroids: Mbh oc a x , w ith x ~ 4.0-4.5 ( Tremai ne et al.l 
2002; Ferrarese & Ford 2005). M BH also increases nearly lin- 
early with galaxy spheroid mass (e.g., Haring & Rix 2004). It 
is now becoming clear that many galaxies have nuclear star 
clusters (NCs) with masses similarly connected to the host 
properties. 

A recent HST/ACS survey of 100 early-type galaxies in 
Virgo has found that ~ 70%-80% of systems with -20.5 < 
Mg < -15 contain dense nuclear components that are resolved 
by HST into star clusters with luminosities L ~ 10 5 L Q - 
5 x 10 7 L©, half-light radii of about 4 pc (possibly with a 
size-luminosity relation among the brighter nuclei), and col- 
ors sugges ting ages of ~2- 10 Gyr and metallicities [Fe/H] ss 
-0.5 ± 1 ICote et a lj2006l). Most spiral bulges and bulgeless 
disk galaxies also contain nuclear clusters with rather simi- 
lar properties dPhi llins et alll996UC arpllo. Stiavelli. & Mackl 
[T99alB5ker et al 120021: 1 Walcher et alJ2005l) . 

All known NCs have masses less than a few 1O 8 M , while 
most measured SMBHs have Mbh > 10 8 M Q . The apparent 
paucity of low-mass SMBHs is at least partly a selection ef- 
fect, bu^the_upperhmit on the nuclear clusters may well be 
real. lCote et alJfeOOfl) find that NC luminosity increases with 
spheroid luminosity, such that Mnc >2x 10 8 M q is expected 
for galaxies with Mb < -20.5. But although they have HST 
surface ph otometry of the cores of all such galaxies in Virgo, 
ICoj^y^ find no evide nce for nucleation in any of them. 

iFerrarese et alJd2Q06h have obtained long-slit spectra for 29 
of the nucleated ellipticals in Virgo. They find that the NC 
masses correlate well with the galaxies' velocity dispersions 
averaged over an effective radius (« 1-2 kpc): Mnc oc a x 
with x = 4.3 ± 0.6, essentially the same as for the SMBH rela- 
tion. However, they also find an offset between the cluster and 
black-hole scalings. Fitti ng power-law M-a relations with 
x = 4 to the Ferrarese et al. NC data and to SMBH data from 
the literature yields (L. Ferrarese, private communication) 



logM NC , 8 = (1.25 ±0.55) + 4 log CT200 
logM BH ,8 = (0.25 ±0.33) + 41ogcT 20 o 



The nuclei of spheroids with velocity dispersion less than this 
are dominated by stellar clusters, with the mass of any SMBH 
that might be present expected to be ~ 10 times smaller. More 
massive galaxies apparently always contain nuclear SMBHs 
b ut, so far as is known, n ot NCs. 

IWehner & Harris! (|2006j use photometric data in the lit- 
erature for about 40 dwarf-elliptical nuclei to show that 
Mnc increases almost linearly with galaxy spheroid mass. 
IFerrarese et alJ J2006I) also find this for their nucleated Virgo 
galaxies with dynamical mass estimates. Moreover, both 
studies conclude that the MNc-M sp h er and MBH-Mpher rela- 
tions meet almost seamlessly at a mass scale ~ 10 8 M Q , i.e., 
there is no large offset as in the M-a relations. These authors 
therefore refer to nuclear clusters and supermassive black 
holes together as "central massive objects," or CMOs. We 
adopt this term here. 

A derivation of the Mbh-c scalin g has been 
given bv iKinl J20031 Eoo3l (see also iFabi^l IT9991 
Murra y. Ouataert. & Thompsonl 120051 iBegelman & Nathl 
2005). He considers super-Eddington accretion on to a seed 
SMBH at the center of an isothermal dark-matter halo. Ac- 
cretion feedback produces a momentum-driven superbubble 
that sweeps ambient gas into a thin shell, which expands into 
the galaxy. Eventually the shock cooling time becomes so 
long that the shell becomes energy-driven and accelerates to 
escape the galaxy. This truncates accretion and freezes in a 
relation of the form Mbh oc a 4 that, with no free parameters, 
matches the observed one remarkably well. A roughly linear 
relation between Mbh and M sp h e r is also established as part of 
this process. 

In this Letter, we examine the possibility that the mass of 
a central star cluster in a protogalaxy might be similarly self- 
regulated, by feedback from stellar winds and supernovae. 1 

2. THE Mcmo — o~ RELATION 



(1) 



where M 8 = M/10 8 M Q and cr 20 o = cr/(200 kms '). A limit 
of M NC < 2 x 10 8 M Q thus corresponds to a < 120 km s" 1 . 



The argument of King ( 2003, 2005) for the M BH -a relation 
has gas in a protogalaxy flowing in to a low-mass, seed black 
hole at super-Eddington rate. This grows the mass of the hole, 
but also drives an intense outflow with momentum flux given 

1 After this paper was submitted, a preprint appeared by 
Li. Haiman. & Mac Low 1 2006), which uses numerical simulations to 
study the development of an apparently common McMO-^snher relation for 
both NCs and SMBHs. Unlike that work, in our analytical model we consider 
stellar feedback explicitly to explain this result, and moreover we attempt to 
understand the offset between the M-a scalings for NCs vs. SMBHs. 
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by the Eddington luminosity: Mv w ~ L^m/c, inde pendent of 
the actual supercritical accretion rate (King & Pounds 2003). 
Here v w is the outflow velocity and Lead = 4-itGMbhc / k, with 
k = 0.398 cm 2 g" 1 the electron scattering opacity. 

King (2003) shows that this outflow is initially momentum- 
conserving, as the shocked gas cools efficiently. The am- 
bient medium is swept up into a thin supershell, which is 
driven outwards by the ram pressure of the SMBH wind: 
p w vl = Mv w /(4nR 2 ) = GM bh /(kR 2 ) at radius R. The dark- 
matter halo is assumed to be an isothermal sphere and the 
ambient gas fraction spatially constant, so that p dm b(R) = 
fgO 2 /(2irGR 2 ). Realistically, there could be a gradient in f g 
due to a concentration of cool gas towards the center of the 
galaxy, but here we take f g to be everywhere e qual to its av- 
erage over the entire halo: f g = Qt,/n m = 0.16 llSpergel et al.l 
2003). If gravity is ignored, the supershell accelerates once 
M B h has grown to the point that p u .v 2 > p am bcr 2 . Any am- 
bient gas outside it is then driven out of the galaxy, stop- 
ping the growth of the SMBH. This happens when Mbh = 
f gK a 4 /(27rG 2 ). 

Including the gravity of dark matt er inside the superbubble 
alters this result by a factor of two (King 2005). The initial 
dynamical expansion of the shell then stalls at a radius /? sta n ~ 
(1 — Mbh/ Merit) -1 , with 

M cvil = f gK a 4 /(nG 2 ) , (2) 

where the gravity balances ram pressure. So long as Mbh <C 
M cr it, this happens well inside the galaxy. More gas can then 
filter through to the nucleus, feeding the hole and causing 
the shell to re-expand (on a Salpeter timescale) to a larger 
i? sta u appropriate to the new Mbh- As Mbh approaches M c ,;t, 
however, the stall radius becomes very large, and before the 
shell can actually reach it the gas cooling time becomes longer 
than the crossing time of the bubble. The shell then enters an 
energy-conserving snowplow phase and accelerates to escape 
the galaxy, leaving Mbh — M cl -it- As emphasized above, it is 
noteworthy that M c ,;t contains no free parameter. 

Our main point is that the above argument is qualitatively 
unchanged if the CMO is not an SMBH but instead a very 
young star cluster where massive stars are still present. Then 
stellar winds and supernovae drive a superwind from the nu- 
cleus with a momentum flux that is much less than L^/c but 
still directly proportional to it. We can thus treat the two types 
of CMO simultaneously by parametrizing the wind thrust as 

Mv w = A L Eid /c = A (47tGM C mo/k) . (3) 

Here A takes a value ~ 1 in the black hole case, but a value 
<C 1 (related to the mass fraction in massive stars) for a nuclear 
cluster. The limiting mass in equation (2) becomes 

Mcmo = 3.67 x 10 8 Mq A" 1 cr 4 200 (/ g /0.16) , (4) 

and the offset Mbh-c and Mnc-c relations of equation ([3 
follow immediately if A ~ 0.1 for a typical NC. 

3. EFFICIENCY OF STELLAR FEEDBACK 

To evaluate the efficiency A of the massive-star feedback 
from a young nuclear cluster, we rewrite equation Q as 

= Mv u . Mvv 

NC ~ 4ttGM nc /k ~ 4.2 x 10 27 (M NC /M Q ) dyn 



and estimate the separate contributions from supernovae and 
stellar winds. 

First, the combined momentum flux from all supernovae 
is 2/V S n£W(vsnTsn), where N SN « 0.01 1 (M NC / M^) is the 
numbe r of stars with mass >8M Q in a cluster with a Chabriei 
(2003) IMF; £sn = 10 51 erg is the energy released per su- 
per nova; vsn ~ 4000 km s " 1 is the typical ejecta veloc- 
ity dWeiler & Srameld H988T) : and t sn rj 2 x 10 7 yr is the 
main-sequence lifetime of an "av erage" SN progenitor (see 
iLeitherer. Robert. & Drisse"nTl 992). Putting this into equation 
(|3 gives Asn ~ 0.02. 

Second, the line-driven wind from a single hot star pro- 
duces a momentum flux of «(L*/c) on average — somewhat 
less than this if there are few lines to drive the wind, but sev- 
eral times higher for O and Wolf-Ra yet stars in whi ch photons 
are multiply scattered (see, e.g., Lamers & Cas sineihlll999l 
chapter 8). Using the m ain-sequence mass-luminosity rela- 
tion of iTout et alJ {1996) to integrate (L*/c) over all stars 
more massive than 5M Q in the IMF of Chabrier (2003), we 
find that the total momentum flux from stellar winds is about 
1.3 x 10 26 (M nc /Mq) dyn, and thus A winds ~ 0.03. 

Despite its very simple derivation, our final 

Anc = Asn + A„i n ds ~ 0.05 (6) 

is in good agreemen t with the values implied by t he detailed 
calculations of Leitherer, Robert, & Drissen ( 1992) for the to- 
tal momentum deposition in solar-metallicity starbursts. Note 
also that the stellar luminosity corresponding to the limiting 
NC mass in equation © with A = 0.05 is co mparable to that 
derived by Murray, Ouataert, & Thompson (2005) from re- 
lated considerations (see their eq. [18]). 

One caveat here is that, while SN momentum fluxes are 
insensitive to stell ar metallicity, wind momenta are ro ughly 
proportional to Z (Leitherer. Rob ert. & Drissenlll992t) . For 
(Z/Zq) ss 1 /3, typical of NCs in Virgo, this might then im- 
ply a net Anc ~ 0.03. However, this effect could be eas- 
ily balanced by increases in both Asn and A w j n d s if the IMF 
in these dense, central starbursts were slightly "top-heavy," 
as may be the case near the center of the Mil ky Way (e.g., 
iNavakshin & Sunvaevll2005t IStolte et alJl2005T) . We proceed 
assuming the fiducial value for Anc m equation (j6ji. 

With A f« 0.05 for a nuclear cluster and A ~ 1 for a su- 
permassive black hole, equation @ implies an offset of a 
factor of 20 between the two Mcmo-c relations, while ob- 
servationally it is only a factor of 10 (ea. fTl). However, 
population-synthesis models (e .g., lFioc & Rocca-V olmerange 
1 19971 JBruzual & Chariot 2003 ) show that a star cluster more 
than ~ 10 9 yr old will have lost some 40%-50% of its initial 
total mass to stellar winds and supernovae, and to the con- 
version of massive stars into degenerate remnants. Thus, we 
expect the Mnc-c scaling originally to have been more offset 
from the SMBH correlation than it is now, by an additional 
factor of about 2. By the same reasoning, the fact that the two 
McMO-M S pher relations currently appear to have very similar 
normalizations must be something of a coincidence. 

4. THE Mcmo -M SPllER RELATION 

In this feedback-regulated picture of CMO and galaxy for- 
mation, an Mcmo-c relation emerges as the primary corre- 
lation. A relation between Mcmo and spheroid mass follows 
by combining equation Q with details of the cooling of the 
wind from the central object. The basic steps are outlined in 
iKind (120031) . Knowing how the shocked gas cools, we find 
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the cooling timescale f coo i as a function of supershell radius 
R and compare it to the dynamical time ffl ow =/?/v s heii. For 
small radii, f coo i < ffl ow and the outflow is momentum-driven. 
However, when the CMO is at about the critical mass in equa- 
tion (|4j, the bubble is so large that f coo i exceeds ffl ow , the thin 
shell becomes energy-conserving, and the wind can escape 
the galaxy. We use /? coo i to denote the radius at which this 
happens. The detailed fate of the swept-up ambient gas after- 
wards is beyond the scope of this Letter, but in general terms it 
should re-collapse to much smaller radii (since the CMO wind 
that pushed the gas to large R in the first place carries essen- 
tially no angular momentum). It will cool rapidly as it does, 
because f coo i < ffl ow by construction inside R < R co0 \, and ?fl ow 
is of order the free-fall time in the halo. We therefore expect 
most of the ambient gas in the supershell at the point of wind 
blow-out to form a concentrated stellar spheroid, and thus we 
identify the ma ss of the shell at R coo \ with M sp h e i-. 

iKind (120031) shows that for a relativistic wind from an 
SMBH, the swept-up gas cools by Compton scattering, and 
ultimately, 



Mi 



-™ = 1.6 x ID" 3 b*l\c/ v JI\f g /0A6T^M:l^ n . (7) 

Mspher 

Here v w ~ c; b ~ 1 is an outflow collimation parameter, 
and M S ph er is in units of 1O H M . This agrees well with 
the observe d SMBH-to-spheroid mass ratios in giant galax- 
ies (lHaring & Rixl2004l) . 

If the CMO is a star cluster, equation Q no longer applies 
because the wind driving the superbubble is far from relativis- 
tic, and the Compton cooling time for the shocked gas exceeds 
a Hubble time. The cooling in this case is by atomic transi- 
tions. However, the M N c-M sp h er relation still has the basic 

form of equation ([7), i.e., Mcmo °c ^tpLi m a ^ cases - 

To see this, note first that we always have ffl ow = R/ v s h e ii = 
R/V2a for the shell as it escapes. In the SMBH case, the 
Compton cooling time is proportional to the inverse of the 
radiation energy density, which is diluted by the l/R 2 law: 
f coo i oc (L Edd /4irR 2 c)~ l oc R 2 /M B n. In the NC case, the cool- 
ing time has the same dependence because the wind density 
also falls off as 1 /R 2 : p w oc M N c/(47r^ 2 v^), and then f coo i oc 
p~ l oc R 2 /Mnc- Thus, either type of CMO has f coo i = fflow at a 
radius /? coo i oc Mcmo / f , or simply R coo i oc er 3 using equation 
(E). Finally, M sp her oc a 2 R coo \ oc a 5 and h ence Mcmo cx M^,.. 
In detail, the flow time is always (e.g.. lKingl2003l) 



erg cm 3 s" 1 for (Z/Z Q ) = 0.1-1 and T ~ 0.5-5 x 10 6 K 
(v u .~ 200-600 km s" 1 ). 

Finding the radius at which f coo i = fflow leads to M sp h er = 
2 f g (T 2 R coo i/G, and combining with equation (0} gives 



MnC .=2.7x10- 4 A-> 



•^sphei 



z 



-0.48 



u_ 

0.16 



-3/5 



.,4.4 



-1/5 



V w,300-'"spher,ll 

(10) 

Comparing equation @ to equation (llOt with A = 0.05, the 
predicted offset of the original MNc-M sp h er relation from the 
corresponding SMBH relation is only a factor of w 3-4- for a 
wind velocity near 300 km s -1 and slightly subsolar metallic- 
ities. Allowing for the long-term mass loss from the NC, dis- 
cussed at the end of ^3] the normalizations of the present-day 
scalings should agree, rather fortuitiv ely, to within a factor 
of two — just as Ferr arese et alJ J2006I) and IWehner & Harris! 
(2006) infer observationally. 

The contrast between this and the much larger offset sep- 
arating the two Mcmo-c relations is due to the fact that, for 
given <7, the radiative f coo i with an NC at the limiting mass 
of equation is nearly 10 times shorter than the Compton 
cooling time that applies when the CMO is an SMBH (see 
lKingl2 003). Thus, the final R coo \ and M sp h er are larger by this 
amount for a galaxy containing a nuclear star cluster vs. one 
with the same a but a central black hole. 

5. THE FABER-JACKSON RELATION 

For a galaxy with a supermassive black hole in its nucleus, 
equations |7} and (|4) imply a relation be tween the m ass and 
velocity dispersion of the spheroid alone (King 2005): 
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M sphel -(SMBH) = 2.8 x 10 11 M Q b- l {c/v w )- 2 (f g /0.\6) 2 a 

(11) 

Perhaps despite appearances, this prediction is consistent 
with the well-known relation between the v elocity dispersion 
and lu minosity of giant ellipticals, L oc er 4 ( iFaber & Jacksonl 
1976). This is because the mass-to-light ratio in the cores of 
these galaxies i ncreases systematically as the ~0.2-0.3 po wer 
of luminosity (van der M arell991t lCappella ri et all2006l) . so 
that the Faber-Jackson relation in fact implies M sp h e i oc a 5 . 
The normalization in equation il It i s also in rema rkably good 
agreement with observation: from Hasegan et al. (2005), 

M spher (obs) ~ 1.93 X 1O U M a" (12) 



fflo W = 6.6 x 10 6 yr fl kpc a 200 A^M" ' (fg/0- 16) 



,1/2 



(8) 



for the shell radius R in units of kpc and M« = Mcmo / 1 s M Q . 
Again, A = 1 for an SMBH and A w 0.05 for an NC. In the 
latter case, we further find for the radiative cooling time, 



'cool 



1 .4 x 10 4 yr R 2 kpc v^ 00 A" 1 Mg 1 (Z/Z )" 



0.6 



(9) 



where v„, 3oo is the speed of the cluster superwind in units 
of 300 km s _1 . This follows from the definition f coo i = 
fimnkT /(p w Aff), with p, ~ 0.6; the wind density p w given 
by the continuity equation, 4TrR 2 pv 2 , = Mv w ; the shock 
temperature kT / pm H = (3/16 )v 2 ,; and Ay the nor malized 
cooling function calculated by ISutherland & DopTtal ( 119931) . 
This last is approximately A N ~ 3.55 x 10" 18 (Z/Z Q ) - 6 r-°- 75 



for galaxies with a > 100 km s" 1 — which, as we discussed in 
^ are the ones that contain SMBHs rather than NCs. 

For smaller galaxies with nuclei dominated by star clusters, 
fcooi = fflow gives the cooling radius 



#cooi(NC) = 880 kpc (Z/Z ) uo C^/O.16)v^oo^2OO • (13) 

As we mentioned above, this is nearly 10 times larger than 
the equivalent scale in the SMBH-dominated case, and the 
spheroid mass is consequently larger than in equation (II 11 1: 

M spher (NC) = 2.6 x 1O 12 M (Z/Z ) a6 (/ g /0. 16) 2 v- 5 ^ 4 )0 . 

(14) 

In a plot of a vs. M sp h er , we would therefore expect low- 



mass, nucleated galaxies to define a locus with a oc 



0.2 
spher' 
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parallel to the non-nucleated galaxies with SMBHs but falling 
below them by a factor of w 10~ 2 . Then, although our 
model does not predict a value for the final effective radius 
of the stellar spheroid in terms of R coo \, the virial theorem re- 
quires any such scale to depend on the spheroid mass roughly 
as R e ff oc M^ er , with nucleated galaxies lying above non- 
nucleated ones by a factor of about 2.5. 

6. DISCUSSION 

Current data probably do not rule out the idea that nuclei 
of galaxies with u < 120 km s" 1 could harbor both star clus- 
ters and SMBHs some ~ 10 times less massive, and thus it 
will be important to ask how such galaxies might choose be- 
tween SMBH and NC feedback channels in regulating their 
formation. For now, a possibly more straightforward question 
is why there are no nuclear clusters in larger galaxies with 
a > 120 km s" 1 (M sphel - > 2-3 x 1O 1O M ), which apparently 
all contain black holes. 

For our basic scenario to be self-consistent, we evidently 
require that the dynamical time of the superbubble always be 
shorter than that of the entire halo, until the point of blow-out 
when ffl ow = f coo i. Equivalently, the radius of the shell must 
always be less than the halo virial radius, and in particular 
^cooi ^ ^vii is required for the growth of a CMO (whether an 
SMBH or an N C) to be self-regu l ated a s we envision. From 
the relations in Brvan & Norman ( 1998), if fl m .o = 0.3, Qa = 



0.7, and H = 70 km s" 1 Mpc" 1 , then R Vlr w 540 kpc a 2 oo (1 + 
z) (accurate to better than 10% for z < 2). Combining this 
with equation Jl 3i . our picture can work for NCs forming at 
redshift znc only in halos with a < 160 km s -1 (1 + Znc)~° 55 ■ 
If nuclear clusters are typically ^5 Gyr old, znc — 0.5 and 
this upper limit becomes a < 130 km s" 1 . In halos with 
velocity dispersion higher than this, the superbubble blown 
by an NC reaches the halo virial radius before it becomes 
energy-conserving and can accelerate to escape. It is pre- 
sumably then held there, or even driven into collapse, by 
the infall of material from beyond R ml . The growth of the 
central star "cluster" is never choked off, and it ultimately 
becomes indistinguishable from the galaxy spheroid itself. 
When the CMO is a b lack hole, th e much longer Compton- 
cooling time found in lKind (|2003) implies that R coo \ < R w i r 
for all a < 500km s _I (1 + Zbh)~ 0,55 , so it is only in very mas- 
sive systems indeed that the self-regulated Mbh-c relation of 
equation (|2j breaks down. 
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Note added in proof. — J. Rossa et al. (AJ, 132, 1074 [2006]) have recently analyzed a sample of nuclear clusters in 40 spiral 
galaxies, finding a relation between Mnc and the luminosity of the galaxies' bulge components, which has the same slope as the 
relation between SMBH mass and bulge luminosity with a comp arabl e intercept. The situation for NCs in late-type spheroids 
thus appears very similar to that found by Wehner & Harris! (|2006) and F errarese et alJ J2006I) to hold in early-type galaxies. 



